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O ■ ABSTRACT 

to: 

, Context. Short time-scale radio variations of compact extragalactic radio sources, known as IntraDay Variability (IDV), can be explained in at 
i ^ | ' least some sources by a source-extrinsic effect, in which the variations are interpreted as scintillation of radio waves caused by the turbulent 
\ interstellar medium of the Milky Way. One of the most convincing observational arguments in favour of propagation-induced variability is the 
t-H ■ so called "annual modulation" of the characteristic variability time-scale, which is due to the orbital motion of the Earth. So far there are only 
^ ' two sources known, which show such a well-defined seasonal cycle, a few more sources with fewer data can be regarded as possible candidates 
for this effect. However, source-intrinsic effects, such as structural variations, can also cause the observed changes of the variability time-scale. 
Data for the new, recently discovered and highly variable IDV source J 1 128+5925 are presented. 

Aims. We study the frequency and time dependence of the IDV in this compact quasar. We measure the characteristic variability time-scale of 
the IDV throughout the year, and analyze whether the observed changes in the variability time-scale are consistent with annual modulation. 
■ Assuming a radio wave propagation effect as origin, we are able to constrain some physical properties (such as distance, scattering-strength 

and possible anisotropy) of the "plasma" screen, which may cause the scintillation. 
. Methods. We monitored the flux density variability of J 1 128+5925 with dense time sampling between 2.7 and 10.45 GHz. We observed with 
O ■ the 100 meter Effelsberg radio telescope of the Max-Planck-Institut fur Radioastronomie (MPIfR) at 2.70 GHz, 4.85 GHz, and 10.45 GHz, and 
J> with the 25 meter Urumqi radio telescope (China) at 4.85 GHz. From ten observing sessions, each of several days in duration and performed 
• ^ ■ between 2004 - 2006, we determine the variability characteristics and time-scales and investigate these in view of possible scintillation and 
p\ ' annual modulation. 

Results. The observed pronounced changes of the variability time-scale of J 1128+5925 are modelled with an anisotropic annual modulation 
model. The observed frequency dependence of the variation is in good agreement with the prediction from interstellar scintillation. Adopting a 
simple model for the annual modulation model and using also the frequency dependence of the IDV, we derive a lower limit to the distance of 
the scattering screen and an upper limit to the scintillating source size. The latter is found to be consistent with the measured core size from 
Very Long Baseline Interferometry (VLB I). 
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1. Introduction 

The so-called IntraDay Variability (IDV) of compact flat- 
spectrum radio sourc es (quasars and blaz ars) was discovered 
in the mid-eighties by IWitzel etaL dl986l) and iHeeschen et al. 
dl987h . If interpreted as being source intrinsic, the short 
variability time-scale (of 1-2 days) would imply - through 
causality and light-travel time arguments - extremely small 
sizes of the emitting regions (yuarcsecond-scale), leading to 
very large apparent brightness temperatures (in the range 
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of 10 16 Kto 10 21 K). This is far in excess of the inverse- 
Compton limit of ~ 10 12 K dKellermann & Pauliny-Tothl 1969b 
and would require with the assumption of a spherical ge- 
ometry excessively large Doppler boosting factors (6 >> 50). 
In somewhat special, non-spherical sourc e geometries (su ch 
as cylindrical shock-in -jet models e.g. lOian et all Il996albl 
199ll:ISpada etdll999h ." the requirements for the Doppler fac- 



tor would be less extreme. Alte rnatively, coherent processes 
and c ollective plasma emission dBenfordlll992l iLesch & Pohll 



19921), or a short t ime intrinsic violation of the inverse Compton 



limit ( Slvshlll992 ) are processes used to explain IDV. 
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While so far no direct other observational evidence exists 
for the assumption of such extremely large Doppler factors, the 
interpretation of IDV via source extrinsic effects appears to be 
very successful, and in particular explains the very fast (hourly) 
variations seen in some extreme ly rapid IDV sources (e.g.: 
Dennett- Thorpe & de BruyrfeOOOh . [n this source-extrinsic in- 
terpretation (e.g.: iRickettl [l990h . the observed IDV is com- 
pletely caused by interstellar scintillation (ISS) of radio waves 
passing through the turbulent plasma of the Milky Way (prop- 
agation effect). However, the emission and source structure of 
compact extragalactic radio sources (quasars, blazars) is known 
to vary also intrinsically on almost every time scale, from min- 
utes to years. Thus it is very difficult to discriminate between a 
source intrinsic and a source extrinsic origin of IDV at radio- 
bands. 

One of the most convincing arguments in favour of a 
source extrinsic explanation of IDV is the so called annual 
modulation of the variability time-scale. An annual cycle in 
the observed variability time-scale is caused by the varia- 
tions of the relative velocity vector between the scattering 
screen and the v elocity of the Earth as the Earth orbits 
around the Sun ( De nnett-Thorpe & de Bruynl 120031) . Such 
seasonal cycles are seen in t wo IDV sources: J 1819+3845 
(jDennett-Thoroe & de Bruynl l2003h 
( BignalletaLl |2003|). In a few of 
such as PKS 1519 - 273 dJauncev et alJ 120031) B 0917+624 
dRickett et al] l2Q0lk Ijauncev & Macauartl l200lb. PKS 0405- 
385 dKedziora-Chudczerll2006h and S5~0954+658 (Fuhrmann 
et al. in prep., 2007), the observed variability timescales do 
not show such a clear effect, and a possible seasonal pattern 
is either not present or is smeared out. For two sources, 
B 0917+624 dFuhrmann et alJ l2002h and PKS 0405-385 



scales of hours or less than an hour. In contrast, the variability 



time- scales of classical type II IDV sources dHeeschen et al 
19871) are much longer, typically two days, and so time-delay 



and 
other 



PKS 1257-326 
IDV sources, 



so called episodic IDV (e.g. Kedziora-Chudczer 20061) is 
observed, where previously observed pronounced IDV, ei- 
ther temporarily disappears (PKS 0405-385), or even ceases 
(B 0917+624). This makes it difficult to prove beyond doubt 
the existence of any annual modulation pattern. We note 
that such episodic IDV can also be due to changes of the 
source structure (i.e. expansion of previously scintillating 
structure components) or due to changes in the properties of 
the scattering plasma. 

The other observational argument for an extrinsic explana- 
tion of IDV is the "time-delay measurements", where a time 
delay is measured between the IDV pattern arrival times at 
two telescopes. Such time delay of the variability pattern was 
observed in PKS 0405-385 between the Australia Telescope 
Compact Array (AT CA) and the Very Long Array (VLA) 
Jauncev et ail d200lb . in J 1819+3845 between the Westerbork 
Synthesis Radio Telescope (W SRT, Netherlands) and the VLA 
dDennett-Thorpe & de Bruynll2002l). and in PKS 12 57-326 be- 
tween the ATCA and the VLA dBignall et alJl2006l) . These de- 
lay measurements, when combined with the observations of 
the annual modulation cycle yield constraints to the velocity, 
geometry and the distance of the plasma causing the scintilla- 
tion. The time delay measurements are feasible only for rapidly 
varying IDV sources, where the variability time scale is signifi- 
cantly shorter than the duration of the observing interval. These 
so called fast scintillators have characteristic variability time- 



measurements cannot be conducted for these sources. 

To establish convincingly an annual modulation pattern, the 
source has to be observed regularly during the course of the 
year and for several years. To facilitate the measurement of a 
characteristic variability time-scale, the duration of an individ- 
ual IDV session must be long enough, so that several "scintilla- 
tion" events ("scintles") can be observed. Within the schedul- 
ing constraints of large observatories, the regular observations 
of fast scintillators, which vary on time-scales of a few hours 
or less, is easier to perform than for the slower, classical type II 
IDV sources. For these sources throughout the year regularly 
performed IDV observations (each lasting at least 3-4 days) 
cannot be performed with an oversubscribed telescopes, like 
the MPIfR 100 meter Effelsberg radio telescope. The search 
for annual modulation in classical type II sources has to in- 
voke other radio telescopes, which can provide more observing 
time. 

In this paper, we present ten epochs of IDV observations 
for a recently discovered and rapidly variable radio source 
GB6 J 1 128+5925 (henceforth referred to as J 1 128+5925): six 
epochs measured with the MPIfR 100 meter Effelsberg ra- 
dio telescope at wavelengths 2.8 cm, 6 cm and 11 cm, and six 
epochs measured with the 25 meter radio telescope in Urumqi 
(China) at 6 cm wavelength. These observations were aimed 
to study the IDV characteristics and to measure the IDV time- 
scale throughout the year. 

J 1128+5925 is a radio qu asar at a redshift z = 1.795 



( So wards-Emmerd et alJ 120051) . the 2000.0 epoch equatorial 
coordinates of J 1 128+59 25 are a 2 ooo.o = 1 l h 28 m 13.3406 s and 



520oo.o = +59°25' 14.798" dFev et alJl2004 . 

The radio spectrum of J 1128+5925 was observed at sev- 
eral dates and is shown in Fig.Q] The spectrum peaks between 
8 GHz and 10 GHz. Between 0.365 GHz and 8.4 GHz, a spec- 
tral index based on n on-simultaneous data from the literature 
dPatnaik et al.ldl992l):lGregorv & Condon! dl 99 il):lBecker et alJ 
il99ll) : IWhite & Becker ! dl992h and iDouglas et al] dl996lV ) of 
q 'o 365 I ghz = 0-18 ± 0.07Ljis measured. Quasi-simultaneous ob- 
servations performed with the Effelsberg radio telescope and 
the Plateau de Bure interferometer of the Insititut de Radio 
Astronomie Millimetrique (IRAM) in 2001 yield a flat high fre- 
quency spectrum, with spectral index i^iosghz ~ _ 0.07 + 0.01. 
Quasi-simultaneous measurements performed in 2006 show 
some spectral variability when compared to earlier observa- 
tions, now yielding o^ mghz = ^-07 + 0.04 below the turnover 
frequency, and a^ 2 ^ = -0.22+0.04 above this turnover. The 
spectral shape and variability is consistent with an inhomoge- 
neous synchrotron self-absorbed compact radio source, which 
undergoes some long-term variability. 

In the following we describe the organization of the paper. 
In Sect. 12 we summarize the observational details, in Sect. [3] 
we describe the basic methods used in the variability analysis 
of the data, in Sect. H] we present the results and discuss them 



The spectral index is defined as: S ~ v + 
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• Quasi-simultaneous at 29. 1 1 .2006 
Quasi-simultaneous at 27.07.2001 
Non-simultaneous in 1990s 
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Fig. 1. The radio spectrum of J 1128+5925. G rey diamonds 



summ ari ze flux density m ea surem en ts from |P atnaik et al 



dl992h: iGreeorv & Condon! dl99lh; iBecker et alj dl99lh: 



White & Becked (119921) and iDouglas et all (119961) . The grey 
dashed line represents the fit to these data points, the corre- 
sponding spectral index «03f^GHz is 0.18 + 0.07. Filled grey 
squares denote for quasi-simultaneous flux-density measure- 
ments with the MPIfR 100 meter Effelsberg radio telescope 
and the IRAM Plateau de Bure interferometer, observed on 
22th July 2001. The solid, grey line represents the fit to 
these data points, the corresponding spectral index a' 



5 GHz 
10.5 GHz 



-0.07+0.015. The filled black-circles show quasi-simultaneous 
flux-density measurements obtained with the Effelsberg tele- 
scope on 29th November 2006. The dashed line is a fit to the 
data up to 8.35 GHz; the corresponding spectral index ^ 64GHz 
is 0.07 ± 0.04. The solid black line is a fit to the data be- 
tween 8.35 GHz and 32 GHz; the corresponding spectral index 
«8 2 3?GHz is "0-22 ±0.04. 



in the framework of the annual modulation scenario. Sect. [5] 
contains a summary and conclusion. 

The following cosmological parameters were used 



throughout this paper: Hq = 71 km s Mpc , £2 



matter — 



0.27 and Q v 



0.73. 



2. Observations and data reduction 

The observations were performed with the MPIfR 100 meter 
Effelsberg radio telescope at centre frequencies 2.70 GHz 
(11 cm), 4.85 GHz (6 cm) and 10.45 GHz (2.8 cm), and with the 
Urumqi radio telescope at 4.85 GHz. The 6 cm receiver, a new 
backend and a new driving program for the Urumqi telescope 
was provided by the MPIfR. A detailed technical description 
of th e 6 cm receiver system and the radio telescope are given in 
e.g.F 



Sun et al. (2007) 



J 1128+5925 has been observed at 10 different epochs, 
which are inhomogenously distributed over a 1.5 year inter- 
val between December 2004 and July 2006. In two epochs 
(December 2005 and April 2006), the Urumqi and Effelsberg 
telescopes observed quasi-simultaneously, to allow for checks 



Table 1. Summary of IDV observations of J 1128+5925. The 
table lists the observing dates (Col. 1), the observing radio 
telescopes (Col. 2, 'E' for Effelsberg, 'U' for Urumqi), the 
observing frequencies (Col. 3), the total observing time (Col. 

4) , and the mean time sampling for the flux measurements of 
J 1128+5925 (Col. 5). Letters shown in the last column (Col. 

5) refer to the corresponding variability curve, shown in Fig. [4] 



Epoch 



R.T. 



V 

[GHz] 



Duration 
[h] 



Sampling 
points/hr 



25-31.12.2004 


E 


4.85 


125 


0.2 


(a) 


13-16.05.2005 


E 


4.85 


67 


0.8 


(b) 


13-16.05.2005 


E 


10.45 


35 


0.8 




14-17.08.2005 


U 


4.85 


67 


0.4 


(c) 


16-19.09.2005 


E 


4.85 


64 


1.4 


(d) 


27-31.12.2005 


U 


4.85 


87 


0.5 


(c) 


29-30.12.2005 


E 


2.70 


29 


1.0 




29-30.12.2005 


E 


4.85 


30 


2.0 


(c) 


29-30.12.2005 


E 


10.45 


30 


1.4 




10-12.02.2006 


E 


2.70 


43 


1.1 




10-12.02.2006 


E 


4.85 


43 


2.0 


(f) 


15-18.03.2006 


U 


4.85 


71 


0.6 


(g) 


28.04-02.05.2006 


E 


2.70 


86 


1.7 




28.04-02.05.2006 


E 


4.85 


86 


1.4 


(h) 


28.04-02.05.2006 


U 


4.85 


94 


1.1 


(h) 


10-13.06.2006 


U 


4.85 


77 


0.9 


(i) 


14-18.07.2006 


U 


4.85 


95 


0.8 


a) 



of the consistency of the flux-density variability seen with the 
two telescopes. The details of the observations are summarized 
in Table [U 

At both telescopes, the flux density measurements were 
performed using the cross-scan technique, in which the tele- 
scope was repeatedly moved over the source position in az- 
imuth (Az) and in elevation (El) direction. Four to eight of 
these slews (sub-scans) in each driving direction formed a scan, 
which, because of the point-like brightness distribution of the 
source, is of Gaussian shape. After subtraction of a linear base- 
line, which removes residual system temperature drifts during 
a sub-scan, a Gaussian curve was fitted to each of the resulting 
slices across the source. The peak amplitude of the Gaussian 
yields a measure of the antenna temperature and of the source 
strength. The on-source integration time from the 4-8 subscans 
was found to be adequate to obtain flux-density measurements 
with a signal-to-noise ratio of SNR > 10 4 at all observing fre- 
quencies, which provides sufficient measurement accuracy. 

The Effelsberg receivers at 6 cm and 2.8 cm are dual-horn 
beam-switch systems, which allow to observe simultaneously 
the source and nearby blank sky (meaning, without the source). 
The immediate subtraction of the two signals removes the at- 
mospheric fluctuatio ns from the source sig nal. For a detailed 



descri ption see e.g. 
(2003). 



Heeschen et al. 019871) and iKraus et al 



In contrast to Effelsberg, the 6 cm at Urumqi measured with 
a single horn only, therefore the fluctuating weather effects 
could not be subtracted. This reduces the measurement accu- 
racy when the weather is not optimal, and in the presence of 
rain or rapidly moving clouds led to data, which could not be 
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Table 2. Primary calibrators and their flux densities at the three 
observing frequencies. 



Name 


S 2.70 GHz [Jy] 


S 4.85 GHz [Jy] 


S 10.45 GHz [Jy] 


3C48 


9.32 


5.48 


2.60 


3C286 


10.56 


7.48 


4.45 


3C295 


12.19 


6.56 


2.62 


NGC 7027 




5.48 





used for the further data analysis (data removed by manual data 
editing). 

In the next step, the amplitudes of the Gaussian fits to the 
individual sub-scans were independently averaged in each driv- 
ing direction (Az, El), and small corrections for the residual 
telescope pointing errors were applied. In the last step the so 
corrected amplitudes from both driving directions were aver- 
aged together, yielding a single antenna temperature measure- 
ment for each source and scan. 

The next step of the data reduction is to correct for the sys- 
tematic elevation dependent gain variations. Elevation depen- 
dent effects mainly originate from the gravitational deforma- 
tion of the dish, which causes variations of the telescope gain 
with changing elevation. A polynomial correction was deter- 
mined using the combined gain-elevation curves obtained from 
the non-variable calibrator sources. 

After this, the time-dependent gain variations must be cor- 
rected. These time-dependent effects are mainly caused by 
the changing weather conditions during the observations, by 
long-term gain fluctuations of the receivers and from thermal 
gradients in the telescope. These effects can be corrected us- 
ing a gain-time transfer function derived from measurements 
of several secondary calibrator sources (e.g. B 0836+710 and 
B 095 1 +699), which are observed with the same duty cycle as 
the target sources. Since these sources do not vary on short 
(IDV) timescales, their measured amplitude variations can be 
used to determine the time dependent gain variation of the 
whole system (telescope + receiver + atmosphere) with high 
(< 0.8 %) accuracy. The observation of more than one of these 
secondary calibrators ensures the correctness of the assumption 
of non-variability of each secondary calibrator. 

In the last step, the measured flux densities were then tied 
to the absolute flux-density scale, which was determined from 
repeated observations of the primary calibrators e.g. 3C48, 
3 C 286, 3C 295 and N GC 7027 and using the flux-density scale 
of lBaars et aD (1 19771) and btt et al] d 19941) . In Table |2] we sum- 
marize the used flux densities for each primary calibrator. 

The residual scatter in the measured flux densities of pri- 
mary and secondary calibrator provides a very good and also 
conservative estimate of the overall calibration accuracy. We 
note that this is mainly limited by the systematic time varia- 
tions of the overall system response, and not by the individ- 
ual flux measurement, which is of high SNR. These time de- 
pendent gain variations are monitored through the observed 
variations of the primary and secondary calibrator fluxe den- 
sities. After correction of all afore-mentioned gain effects, the 
residual variability of the calibrators yields a very reliable es- 
timate of the measurement accuracy, and furthermore does not 



depend on a-priori assumptions on receiver stability and tele- 
scope gain. We quantify the overall measurement accuracy us- 
ing the variability index mo =< (100 • crj < S >) > of the 
secondary calibrators, where mo is the weighted average ra- 
tio of the standard deviations to the mean flux density for all 
calibrators. Thus, the uncertainty of an individual flux density 
measurement is be given by mo • S . The values for mo for each 
experiment are given in Col. 5 of Table [3] Typically, one ob- 
tains for Effelsberg m ~ 0.4% at 2.70 GHz, m ~ 0.5% at 
4.85 GHz, and m ~ 0.6 % at 10.45 GHz, mainly reflecting the 
fact that the atmosphere limits the measurement accuracy to- 
wards higher observing frequencies. 

At Urumqi one obtains at 4.85 GHz typically mo = 0.5 - 
0.7 %. A somewhat larger value of mo = 1.2 % is observed in 
December 2005, which is due to bad weather conditions (rain 
and snow). 

2.1. Comparing the Effelsberg and the Urumqi data 

The direct comparison of the variability curves obtained with 
the 100 m Effelsberg and 25 m Urumqi telescope in two quasi- 
simultaneous observing sessions (Dec. 2005, April 2006, see 
Fig. |2] and Tab. [3}, shows a very good agreement of the vari- 
ability curves measured at both telescopes and demonstrates 
that the 25 m Urumqi telescope is well suited for IDV studies. 
The small differences seen in the rms-scatter of both data trains 
are likely to result from some fundamental differences between 
the two observatories: 

- The 100 m Effelsberg dish uses the homology principle, the 
25 m Urumqi telescope does not. Consequently, the gravi- 
tational deformation as a function of elevation is different 
for both dishes, resulting in a lower aperture efficiency and 
steeper gain-elevation curves of the Urumqi telescope. 

- While for the 100 m Effelsberg telescope, the focus posi- 
tion is automatically adjusted at each elevation and also al- 
lows manual displacement and refocusing at any time (day, 
night), the focus in the Urumqi telescope is at a fixed po- 
sition, not allowing for additional focus corrections. This 
leads to stronger variations of the overall time dependent 
gain. 

- At 4.85 GHz, Effelsberg uses a dual-horn beam-switch re- 
ceiver, while at Urumqi a single horn total power receiver 
is used. The Urumqi data are therefore more affected by the 
weather than the Effelsberg data. At Urumqi, this results in 
a less stable signal, which is probably the main reason for 
the larger scatter seen in these data. 

For the 6 cm observations performed at Urumqi, we used 
the following secondary calibrators, each of which with a 
flux density comparable to that of our program sources: 
B 0836+710 (5 4 .85GHz = 2.1 Jy), B 0951+699 (S 4 .8 5 ghz = 
3.3 Jy), B 1203+645 (S 4 .85GHz = l-2Jy), B 1128+455 
(S 4.85 ghz = 0.7 Jy), B 1311+678 (S 4 .85gh z = 0.9 Jy). We find 
it remarkable that despite a 16 times smaller collecting area 
of the Urumqi telescope, the achieved measurement accuracy 
(m = 0.5 - 1.2 % at Urumqi, m = 0.4 - 0.5 % at Effelsberg, 
see Tab. [3} is only ~ 2 - 3 times worse than at Effelsberg. 
This demonstrates that the IDV measurements performed with 
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December 2005 




1.5 2 2.5 

Days since J.D. 2453853.935 



Fig. 2. 4.85 GHz radio variability curves of J 1128+5925 ob- 
tained in December 2005 (top) and April 2006 (bottom). Filled 
grey circles denote observations with the Effelsberg telescope, 
black open circles denote observations carried out with the 
Urumqi telescope. 



Urumqi are not limited by the size of the dish and the signal- 
to-noise ratio, but by the larger fractional errors which are most 
likely the result of the afore-mentioned three effects. 

3. Data analysis 

For the variability analysis, we follow earlier work and use a 
parameterization already introduced by Heesch en et alj (119871) 
and described in detail in e.g. Ouirrenbach et al.h 2000h . 

To decide whether a source can be regarded as variable, we 
performed a % 2 - test, which tests the hypothesis whether the 
light curve can be modeled by a constant function. Sources for 
which the probability of this hypothesis is less than 0.1% are 
considered to be variable. 

The amount of variability in a light-curve is described by 
the modulation index and the variability amplitude. The modu- 
lation index is defined by m = 100 • where (S) is the time 
average of the measured flux density, and <x is the standard de- 
viation of the mean flux density. (The mean modulation index 
of the calibrators, mo describes the overall calibration accuracy, 
see Sect. 12). 

In estimating the errors of the modulation index of the vari- 
able source, the statistic al effects o f a stochastic process have to 
be taken into account dCim6ll2003l) . The errors originate from 
the limited observing time and from the time sampling. If the 
observing time is much shorter than the characteristic variabil- 
ity time-scale (no clear maximum followed by a minimum can 
be observed in the light-curve), the variability index m is not 
well defined. The contribution to the statistical error, originat- 
ing from the length of the observing interval, is proportional to 
V^var/^obsi where f var is the characteristic variability time-scale 
and r o b s is the total observing time (B. Rickett, priv. comm.). 
In addition, variability may also exist at times below the mean 
sampling interval of the measurements. Under-sampled vari- 



ability also leads to an underestimate of the variability index m. 
Assuming that significant variations may exist between three 
data-points, the error is gi ven as: 3 ■ ( t s )/T„y, s , where (t s ) is 
the mean sampling interval dCimoll2003l) . Both of these effects 
were taken into account for the calculation of the uncertainty 
of the modulation index m, which was found to range between 
10 % and 30 % of the given value of m. 

Following Heeschen et al. (1987) we define a 3 cr- variablity 
amplitude Y through the noise-bias corrected 3 cr-value of the 
modulation index m. We define the variability amplitude of a 

given source with variability index m as Y — 3 ^jm 2 - m^, 
where mo is the mean modulation index of all secondary cal- 
ibrators, which describes the statistical measurement accuracy 
during the observation. Y is set to zero for non-variable sources 
(m = mo). 

3.1. Determination of the characteristic variability 
time-scales 

There are many different definitions for the characteristic 
variability time-scale of IDV used in the literature. E.g. 



Simonetti et alJ (I1985I) use the structure function of the ob- 
serve d time series t o defi ne a characteristic variability time- 
scale. [fadsetteTaj] dHH) and iRickett et alj d2002l) define the 
variability time-scale by t he half-width at half maximu m of the 
autocorrelation function. Macquart & Jaunceyl d2002h use the 
half-width at l/e of the maximum of the autocorrelation func- 
tion (this definition of the so called decorrelation ti me-scale 
is ad opted from pulsar scintillation studies, see e.g. ICordes 
19861) . Since IDV light-curves appear often quasi-periodic, 
resembling some sort of sine-wave pattern, the variability 
time-s cales can also be estimated f rom the mean peak-to 



peak dKedziora-Chudczer et al. 1997 ) or peak-to-trough time 



dJauncev & MacquarfcOQll) . 



To obtain reliable estimates of the variability time-scales 
present in our data, we derived the time-scales using 3 different 
methods: 

- We calculated from the light-curves the average peak-to- 
trough time (t\ c , Col. 12 in Table [3). The errors are cal- 
culated from the internal scatter of the different peak-to- 
trough times found in one data train. In some light curves, 
the peaks and troughs were not well-defined. In these cases 
we do not give an error estimate, the time-scales are given 
in parenthesis in the table. 

- We calculated the structure function (SF) for each time 
series and derived the "saturation" time-scale from it. 
(tge, Col. 13 in Table [3). The SF is defined as SF(r) = 
<(S (t) - S (t + t)) 2 >, where S (f) is the flux density time- 
series, t is the time-lag and () denote time averaging. Above 
the noise level, the SF rises monotonically and is described 
by a power law. At large time-lags, the SF reaches its maxi- 
mum at a "saturation" level, which is proportional to 2 • m 2 . 
The intersection of the power-law fit with this "satura- 
tion" level defines th e characteristic variability time-scale t 
dBeckert et al ] |2002|) . The measurement error for this time- 
scale is derived from the formal errors of the power law 
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fit and the fit of the "saturation" level dFuhrmannl 12004 
Fuhrmann et al 12007 , in prep.). 



- We also calculated the autocorrelation function (ACF) of 
the tim e series, using the method of Edelson & Krolik 
dl988l) . The ACF is related to the SF via: SF(t) = 
2 (ACF (0) - ACF (t)). Thus, the time-lag derived from the 
point where the SF reaches its maximum corresponds to the 
time-lag at the first minimum of the ACF. The derived time- 
scales (?acf) are given in Col. 14 of Table [3] The error of 
?acf are determined from the scatter in the positions of the 
first minimum of the ACF, when calculated with different 
bin sizes. 

Examples of structure functions (SF) and autocorrelation 
functions (ACF) obtained for J 1 128+5925 are displayed in Fig. 
[3] The functions were derived from the light curve observed at 
4.85 GHz with the Effelsberg radio telescope in 28 April 2006. 
The corresponding light-curve is shown in Fig.|4]in subplot (h). 

4. Results 

In Figure [4] we show the radio variability curves of 
J 1 128+5925 obtained at 4.85 GHz during December 2004 and 
July 2006. For each observing date the normalized flux den- 
sity of J 1128+5925 (top panel) and of a secondary calibrator 
(B 0836+710, panel below) is plotted versus time (in days). In 
the plots, the amplitude scale for J 1 128+5925 and for the sec- 
ondary calibrators are similar, visualizing the prominence of 
the variability seen in J 1128+5925. Here, the peak-to-trough 
amplitude is found to reach 30 %. The light-curves often ap- 
pear quasi-periodic (e.g. in 2005 September subplot (d), 2005 
April subplot (h)) and show some complexity, indicative of the 
presence of more than one variability time-scale. The variabil- 
ity characteristics, time-scales and strength of the variability 
change significantly from epoch to epoch, and sometimes even 
within a 3 - 4 day interval. For example in June 2006 (subplot 
(i)), where a pronounced trough and peak is followed by a much 
shallower variations during the last days of the observations. 

In Table [3] we summarize the results of the observations of 
J 1 128+5925. We list the variablity parameters for all observing 
sessions: Column 1 gives a label, which refers to the variabil- 
ity curve shown in Figure @] col. 2 summarizes the observing 
date (first day, dd/mm/yyyy), col. 3 the observing frequency, 
col. 4, the observatory, col. 5 the mean modulation index of the 
secondary calibrators, column 6 the number of data points ob- 
tained for J 1128+5925, col. 7 the mean (time averaged) flux 
density of J 1128+5925 in this session, col. 8 the standard de- 
viation for the mean flux density, col. 9 the modulation index, 
col. 10 the variability amplitude, and col. 11 the reduced^ 2 . 
For the degrees of freedom considered here, a 99.9 % prob- 
ability for variability is given, when the reduced x 1 exceeds 
1.4 - 2.1. With^- 2 in the range of ~ 4 to ~ 1000, the variabil- 
ity of J 1128+5925 is highly significant at all observing dates 
and frequencies. We note that at 10.45 GHz the variations are 
somewhat lower. A detailed description and discussion of the 
frequency dependence of variability is given in Sect. 14.31 

In columns 12, 13, and 14 of Table [3] we summarize the 
three different estimates of the variability time-scale. The cal- 



culation of the variability time-scale is described in detail in 
Sect. ED 

The typical variability time-scales range between 0.2 - 
1.9 days. For each observing epoch and data set, the three dif- 
ferent methods reveal similar results and show the internal con- 
sistency of the measured variability time-scales. 

Adopting a source-intrinsic interpretations for the variabil- 
ity, we can apply the light travel-time argument and derive - via 
the source size - an apparent brightness temperature. Following 
Mars cher et al. the angular size (in mas) is obtained 



from: 6 S - 3.56 • 10 f var (l +z)6/D h . Here, / val is the vari- 
ability time-scale (in days), z is the source redshift, Z\ the lu- 
minosity distance (in Gpc) and S the Doppler factor. This size 
now can be used to calculate a brightness temperature which, 
in the case of a Gaussian brightness distribution, is given by: 
T B = 1.22 ■ 10 I2 S (1 +z) / (v 2 6> 2 ) K (e.g. lKovalev et~ai1l2006l) . 



Here S is the flux density in Jy, v the observing frequency in 
GHz and 9 the source size in mas. 

Even if one uses the longest measured variability time- 
scale (of 1.7 days as measured in April 2006 at 4.85 GHz), 
one obtains a brightness temperatures far in excess of the 
inverse-Compton limit of 7b ~ 10 19 K. The fastest ob- 
served variations (seen in December 2005 and February 
2006) imply T B ~ 10 2() K. Independent of the use of 
the inverse Compton brightness temperature limit (10 12 K, 
Kellermann & Pauliny-Toth 1969h o r the equipartition bright- 



ness temperature limit (10 H K, IScott & Readheadl 1 1977 



mp 

Readheadl 1 19941) . this leads to Doppler-boosting factors of a 



couple of hundreds to thousand, which are necessary to bring 
the observed brightness temperatures back to these limits. 

4.1. Time dependence of the variability 

In J 1 128+5925, the variability amplitudes and time-scales vary 
at each observing frequency significantly between the different 
epochs (see Fig.|4]and TableO. 

The length and sampling of the observations may influence 
the measurements of these parameters. However, only in some 
cases are the observed differences due to sampling effects. This 
is illustrated in the variability curves of December 2005, which 
despite a different time sampling and and by a factor of two 
different duration of the observation at Effelsberg and Urumqi 
still leads to very comparable variability indices, amplitudes 
and time-scales. 

To search for a possible correlation of the IDV with long- 
term variations of the total flux density, we plot in Fig. [5] 
the mean flux density, the modulation index, the variabil- 
ity time-scale and the spectral index versus observing epoch. 
Measurements at different frequencies performed by different 
telescopes are shown with different symbols. 

The variations of the modulation index and of the variabil- 
ity time-scale do not show a correlation with the total flux den- 
sity. The mean flux densities at 4.85 GHz showed an increase 
by 20 % until February 2006, over a 1 yr time-scale. This in- 
crease is also seen at 10.45 GHz. After this, the 4.85 GHz flux 
density decreased almost to its original value by July 2006. 
Unfortunately, we were unable to find additional total flux den- 
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Fig. 3. Example of a structure function (left) and an autocorrelation function (right) characterizing the variability of J 1 128+5925. 
The data are from Effelsberg observations performed at 4.85 GHz in April 2006. 
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Fig. 5. The mean flux density ((S ) in Jy), the modulation in- 
dex (m in percent), the characteristic variability time-scale (in 
hours) and the spectral index of J 1 128+5925 displayed versus 
observing epoch. Black circles and black solid lines denote ob- 
servations at 4.85 GHz: large open circles denote observations 
with the Urumqi telescope, small filled circles denote observa- 
tions with the Effelsberg telescope. Triangles and dashed line 
denote observations at 10.45 GHz, grey squares with dashed- 
dotted line denote observations at 2.70 GHz. Diamonds denote 
for the spectral index. Open diamonds are for the spectral index 
between 10.45 GHz and 4.85 GHz, filled diamonds are for the 
spectral index between 4.85 GHz and 2.70 GHz. Grey diamond 
is for the spectral index between 10.45 GHz and 2.70 GHz. 



sity monitoring data for J 1128+5925 in the literature, which 
would allow to us to check if and how the flux density varies 
on longer time-scales. Clearly, a future flux density monitoring 
would be helpful. 



The spectral index was calculated for those 4 epochs for 
which we had quasi-simultaneous multi-frequency measure- 
ments. We obtained a spectral indices between 4.85 GHz and 
10.45 GHz in May 2005 and in December 2005, and between 
2.70 GHz and 4.85 GHz in February 2006 and April 2006. We 
also determined the spectral index using all three frequencies in 
December 2005. Within the limited amount of multi-frequency 
data, we do not see pronounced temporal variations in any 
of these spectral indices. The observations, however, indicate 
some change of the spectral slope at 4.85 GHz from steep to 
flat. The data are also consistent with a spectral turnover be- 
tween 8 GHz and 10 GHz, which is already seen in Fig.Q] 

4.2. Annual modulation model 

The variability time-scales seen in J 1128+5925 are signifi- 
cantly different at different epochs (Fig. [4] and Table [3j. At 
4.85 GHz, the changes in the characteristic variability time- 
scale do not correlate with other parameters of the variability 
(see Fig. [5j. However, it is striking that in two observations 
separated by one year (December 2004 and 2005) a very simi- 
lar variability time-scale of 0.3 days is seen. The robustness of 
this short variability timescale is independently confirmed by 
two telescopes (Effelsberg and Urumqi, December 2005). 

In Fig. [6] we plot the variability time-scale versus day of 
the year. Different symbols represent different years of obser- 
vations: stars for 2004, squares for 2005, circles for 2006. From 
Fig. [6] the slowest variations occurs around day 120 (end of 
April to early May, f var w 1.6 day). Although the number of 
data points is sparse, the figure shows another "slow-down" of 
the variability time-scale around day 230 (August, September). 
In the following, we try to explain this behaviour by two sce- 
narios within the framework of the annual modulation model. 

A seasonal cycle of the characteristic variability time-scale 
can be explained in a model of interstellar scintillation (ISS), 
where the scattering material is regarded to be located in a 
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Fig. 4. 4.85 GHz light-curves of J 1128+5925 (upper plots) and of one secondary calibrator (B 0836+7 10,lower plots) from 
December 2004 (panel a) until July 2006 (panel j). Letters in brackets at the top left of each plot refer to Table Q] where the 
mean source flux density and the variability parameters are given. The abscissa displays the observing time [in days], relative to 
the starting date, which is shown in the top up right corner of each plot (see also Table [T)- On the ordinate the normalized flux 
densities are shown, with normalization by the mean source flux density. Filled grey symbols denote Effelsberg observations, 
open black symbols represent Urumqi observations. 



thin plasma screen at some distance from the Earth. The or- 
bital motion of the Earth around the Sun leads to annual ve- 
locity changes, resulting from the vector addition between the 
observer's velocity vector and the velocity vector of the screen. 
Thus annual modulation is a pure geometrical effect. A low 



relative velocity between screen and observer results in a pro- 
longation of the observed variability time-scale. Half a year 
later the two velocity vectors are oriented oppositely. This 
then leads to a large velocity difference and correspondingly 
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Table 3. The variability parameters of J 1 128+5925 measured by the Effelsberg telescope and with the Urumqi telescope between 
2004 and 2006. The letters in Col. 1 refer to the subplots of Fig. [4] Col. 2 shows the starting date of the observation, Col. 3 the 
observing frequency, Col. 4 the observing telescope ('E' for Effelsberg, 'LP for Urumqi). In Col. 5 the modulation index of the 
non-variable secondary calibrator is displayed, in Col. 6 the number of flux-density measurements for J 1128+5925. The mean 
flux density, its standard deviation, the variability amplitude, and the modulation index of J 1 128+5925 are given in columns 7, 
8, 9 and 10, respectively. In Col. 11 the reduced^- 2 (^ 2 ) is shown. In Col. 12, 13 and 14 the measured variability times-scales 
are given. They were derived directly from the light-curves (Col. 12), using structure functions (Col. 13), and autocorrelation 
functions (Col. 14). In some light-curves, the peaks and troughs are not well-defined. In these cases we do not give an error for t\ c 
and the time-scales are given in parenthesis in Col. 12. 
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[GHz] 




[%] 




[Jy] 


[Jy] 


[%] 


[%] 




[day] 


[day] 


[day] 


(a) 


25.12.2004 


4.85 


E 


0.4 


24 


0.570 


0.062 


10.9 


32.6 


649.705 


0.3 ±0.1 




0.3 ±0.1 


(b) 


13.05.2005 


4.85 


E 


0.5 


54 


0.611 


0.013 


2.2 


6.4 


19.932 


(0.8) 


0.8$ 


1.0 ±0.3 




13.05.2005 


10.45 


E 


0.6 


29 


0.730 


0.028 


3.8 


11.2 


11.875 


(1.0) 


0.8 ± 0.2 


1.0 ±0.1 


(c) 


14.08.2005 


4.85 


U 


0.6 


20 


0.682 


0.040 


5.9 


17.5 


4.291 


1.2 ±0.2 


0.8 ±0.5 


1.1 ± 0.2 
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4.85 


U 


1.2 


40 
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0.056 


7.9 


23.4 


25.872 


0.3 ±0.2 


0.4 ±0.2 


0.34 + 0.04 


(e) 


29.12.2005 


4.85 


E 


0.4 


64 
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0.052 


7.2 


21.5 


237.791 


0.3 ±0.1 


0.3 ±0.1 


0.30 + 0.06 
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31 
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6.8 
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E 


1.4 


46 


0.826 


0.027 


3.2 


8.7 


5.340 


0.2 ±0.1 


i +0 - 2 


0.3 ±0.1 


(f) 


10.02.2006 


4.85 


E 


0.4 


96 


0.723 


0.049 


6.8 


20.3 


192.638 


0.2 ±0.1 


0.1 ±0.05 


0.1 ±0.04 
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30.6 
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(1.1 ±0.5)^ 


0.5 ±0.1 
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U 
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40 
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0.038 
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0.5 ±0.1 
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0.5 ±0.1 


(h) 
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72 


0.595 
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a) 
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0.7 
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0.035 


5.8 
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0.6 ±0.3 


6 +03 
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" The time-scale obtained from the SF analysis was not consistent with the light-curve. 

* ? SF are much longer than those derived from the ACF and the light-curve. There is an indication of an additional 
plateau in the SF at a time lag of ss 0.5 day. This agrees well with the other two values of the timescale. 

The time-scales estimations are based upon the variations observed during the first one and half day of the obser- 
vation. During the last two days the flux-density variations were much more reduced (see Fig.f4](i) subplot). 
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Fig. 6. The IDV time-scale of J 1 128+5925 measured at 4.85 GHz plotted versus day of the year and fitted by 6 different annual 
modulation models. Different symbols represent observations performed in different years: stars stands for 2004, squares for 
2005 and circles for 2006. Filled symbols represent observations carried out with the Effelsberg telescope, open symbols for 
observations with the Urumqi telescope. Lines represent time-scales from an interstellar scintillation (ISS) model, assuming 
isotropy (left) and anisotropy (right). The labels refer to six different models, which are summarized in TableHJ 
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to fast variability (e.g. iBondi et alj|l994t iRickett et alj|2001 
Dennett-Thorpe & de Bruvnll2003l) . 

In such a model the free parameters are the components 
of the velocity vector of the screen in the plane perpendicular 
to the line of sight, and the scattering length scale (so). The 
scintillation time-scale as a function of the day of the year (T) 
then is given as: t (T) ~ sq/v (T), where v is the relative velocity 
between the observer and the screen. 

In the left panel of Fig.|6]the three different lines represent 
annual modulation models calculated from three different sets 
of input parameters. The first curve, (a), is the best fit to the 
whole dataset and yields a Xr - 4.6. In an attempt to fit the 
second "slow-down" of the variability time-scale around day 
230, we first excluded the data points at day 120. This gave the 
slowest variability time-scale and thus determines the height 
of the first peak. Thus, we obtained curve (b). The^ of this fit 
was 4.1, however including all the data points yielded^ = 8.3. 
To be able to fit the "slow-down" at day 230 one has to exclude 
also the measurement at day 134. This is illustrated by curve 
(c) The Xr of this reduced dataset was 3.0, however using all 
the data points yielded a. — 13.3. To summarize, the first 
"slow-down" at day 120 is not due to one exceptionally long 
variability time-scale, but is constrained by two measurements, 
one obtained in May 2005 and the other in April 2006. 

The parameters of the different curves are given in Table|4] 
In Col. 1, the letters refer to the corresponding curve displayed 
in Fig. [6] In Col. 2 and 3, the components of the screen velocity 
are displayed in Right Ascension (i>ra) and in Declination (v$). 
In Col. 4, the scintillation length scale is given. Fig. [6] shows 
that this annual modulation model can explain only one slow- 
down of the variability time-scale in one year, but not two. 

The previous model assumed an isotropic scintillation pat- 
tern. However, in a more general scenario, scintillation pat- 
terns could also be anisotropic. This is modelled by introduc- 
ing some ellipticity. Including the anisotropy of the scintillation 
pattern introduces two more model parameters: the angular ra- 
tio of the anisotropy (r) and its position angle (y). The scintil- 
lation time-scale now also depends on the direction of the rela- 
tive velocity vector (between the Earth and the screen) and how 
it "cuts" through th e elliptical scintillation pattern. Following 



Bign all et al.l (2006), the scintillation time-scale is given as the 
function of day of the year (T) by the following expression: 



t(T) = 



s yfr 



^ (r) 2 + ( r 2_ 1)(v(r)xS) 2 



(1) 



where so is the scintillation length scale, v (T) is the relative 
velocity between the scintillation screen and the observer. As 
the Earth orbits around the Sun v varies annually, v denotes the 
absolute value of v. S = (cos y, sin y) is the unit vector, defining 
the orientation of the elliptical scintillation pattern, where the 
angle y is measured from East through North. The major and 
minor axis scale lengths of the scintillation pattern are defined 
through a ma j = so an( l a min = So/ Vr, respectively. We note 
that in this paper the definition of the time-scale differs from 
that used bv lBignall et al.l (12006b . This leads to a factor of ~ 2.5 



In the right panel of Fig.|6]the fits for the anisotropic annual 
models are displayed. The parameters of the three model curves 
are summarized in Table [4] 

The simple isotropical scintillation model can not ex- 
plain both, the "slow-down" of the scintillation time-scale at 
around day 120, and the increase of the variability time-scales 
measured aft e r day 230. The anisotropic scattering model of 
BignalletalJ (2006), however, does reproduce both features 
quite well. Owing to its broader shape and lower modulation, 
this second "slow-down" around day 230 is not as well de- 
termined as the peak around day 120, therefore are the cor- 
responding parameters r and ora are not so well constrained. 

In the case of a point source (i.e. a source smaller than the 
Fresnel angle) the scattering length scale is determined by the 
Fresnel angle and the distance to the scattering screen. The 
Fresnel angle is given as 8p = ~\/A/(2nD), where D is the 
distance to the scattering screen and A is the observing fre- 
quency. Extragalactic radio sources usually cannot be regarded 
as point-like in this context, since 6 S > 6p. Therefore, the scat- 
tering length scale is mainly d etermined by the source size an d 
not by the Fresnel angle (e.g. lGoodmanl ll997; Walker 1 998b . 
Therefore, the distance to the screen can be derived from the 
scintillation length scale with the assumption of an angular 
source size. 

For this, we use the VLBI data obtained for J 1128+5925 
through images and Gaussian modelfits. The source was ob- 
served with VLBI at 5 GHz in t he Caltech-Jodrell Bank Flat- 
Spect rum Survey (CJF survey, Taylor et al. 19961: Ma et al 



differe nce between the scin tillation length scale derived here 
and bv lBignall etail d2006l) . 



19981 and S. Britzen in prep., 2007). J 1128+5925 is a very 
compact radio source, with a possible jet-like feature oriented 
towards the south-west at a position angle of -100°. From cir- 
cular Gaussian modelfits an upper limit to the VLBI core size of 
0.53mas (in 1992.731) and 0.13 mas (in 1994.702) (S. Britzen, 
priv. commm.) is obtained. 

An independent size estimate is obtained from the turn-over 
frequency of the radio spectrum of J 1128+5925 (see Fig. |T}. 
Assuming that the observed spectral turn-over is due to syn- 
chrotron self-absorption and further assuming energy equipar- 
tition between particles and magnetic field, allows to calculate 
the so called equipartition source size (e.g. IScott & Readhead 
Il977b . Adopting 8.35 GHz as the turnover frequency and us- 
ing the flux-density measured at this frequency (0.539 Jy), one 
obtains the equipartition source size via: 6 eq = 0.34D L I/,17 (1 + 
z )9/i7£-7/i7 mas ^ wnere £) L j s the luminosity distance, z the red- 
shift of the source and 6 the Doppler factor. From the CJF sur- 
vey (Britzen et al. 2007, and in prep.) the Doppler factor is esti- 
mated to be 6 ~ 5. With this one obtains an equipartition size of 
» 0.26 mas at 8.35 GHz. We now may assume that the source 
size is inversely proportional to the observing frequency (for 
discussion see, Sect. I4.31 >. This leads to an extrapolated source 
size of 0.45 mas at 4.85 GHz. 

From the definition of the scintillation length scale so = 
6 S D, and using the above derived upper limit for 6 S , we ob- 
tain a lower limit to the screen distance. Taking the source size 
derived from the spectrum of J 1 128+5925 (0.45 mas), the cal- 
culated lower limits on the screen distance range between 10 pc 
and 20 pc. Using the source size derived from the modelfits of 
the CJF survey data (0.13 mas) the lower limits on the screen 
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Table 4. Fit parameters for the annual modulation models shown in Fig. [6] Col. 1 show labels used for the different fit curves. Col. 
2 and 3 give the velocity components of the screen in Right Ascension and Declination direction. Col. 4 shows the scintillation 
length scale, Col. 5 the axial ratio, and Col. 6 the position angle of the scintillation pattern. In Col. 7 we give the reduced ;f 2 for 
the fits. For details of calculation of^- 2 see text. 
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distances range between 36 pc and 60 pc. At these distances 
the Fresnel angular scales at 4.85 GHz indeed range between 
15 yuas and 36 yuas. This confirms our initial assumption of the 
Fresnel scale being smaller then the source size. Therefore the 
source cannot be regarded to be point-like. 

If we assumed that the source size is smaller than or equal 
to the Fresnel scale, than the scintillation length-scale is deter- 
mined by the distance to the scattering screen and the Fresnel 
angular size. However, in that case the assumed source sizes 
would be extremely small (~ 1 /ias) compared to usual extra- 
galactic source sizes. Assuming a circularly symmetric source 
with FWHM corresponding to these Fresnel angle values im- 
plies a lower limit on the brightness temperature of the order of 
10 15 K. 

The modelfits of the CJF survey VLBI data revealed in two 
epochs (1994.702 and in 1996.646) a secondary (jet?) compo- 
nent at a position angle -95° and -104° relative to the core. 
These position angles agree remarkably well with the position 
angle of the anisotropy derived from the anisotropic annual 
modulation model (x -90°). It is therefore possible that the 
annual modulation and the by the model predicted amplitude 
and seasonal times of the two major maxima in Figure [6] are 
related to the VLBI structure and its orientation. At the mo- 
ment it is unclear, if the agreement between both angles is a 
chance coincidence. Moreover the VLBI observations of the 
source were performed ten years before our IDV observations 
and it is unclear if the jet orientation is persistent over such 
a long time. To confirm whether the orientation of the source 
elongation remained constant over time, and whether it is re- 
lated to the anisotropy seen in the annual modulation pattern, 
requires further VLBI observations. 

A seasonal cycles seen in the variability time-scales 
of some other IDV sources (suc h as J 1819+3845 



(iDennett-Thorpe & de Bruviil 120031). PKS 12 57-32 6 



Bignall et al 



2003h 



1519- 273 d.Tauncevetal.1 2003 



B 0917+624 (IRickett et alJ 120011: Uauncey & Macquart 



2001b and S4 0954+658 (L. Fuhrmann priv. com.)) are sim- 
ilarly explained in the framework of the annual modulation 
theory. Assuming anisotro p y in the scintillation p attern , 



Dennett-Thorpe & de Bruynl d2003l) and iBignall et al.1 d2006h 
successfully modeled the annual variations in the characteristic 
variability time-scales in the extremely fast IDV sources, 
J 1819+3845 and in PKS 1257-326. The derived screen dis- 



tances in these sources were less than ~ 10 pc, comparable to 
the lower limits to the screen distance derived here. 

4.3. Frequency dependence of the variability 

J 1128+5925 has an inverted spectrum, with spectral index 
ranging between ~ 0.2 and ~ 0.5 at all epochs (Pig. [5j. 

Fig. [5] also shows, for those epochs with simultaneous 
multi-frequency measurements, a systematic frequency depen- 
dence of the modulation index, decreasing in magnitude from 
2.70 GHz to 10.45 GHz. 

This frequency dependence is shown in Fig. [7] It is con- 
sistent with the ISS theory, which explains the phenomenon of 
IDV by inter s tellar s cintillation in the so called weak regime 



(e.g. iRicketl Il986t IRickett et all 120061) . The variations at 



tributed to refractive scintillation in the stron g regime have 



longer time-scales , usually several days (e.g. IRickett 1986; 



Rickett et al.ll2006l) . The characteristic time-scales of the varia 



tions observed in J 1 128+5925 show a slight increase with de- 
creasing frequency in all of the multi-frequency datasets, which 
is also consistent with the theory of WISS (weak interstellar 
scintillation) (see Fig. [5j. 

In the weak regime, m decreases with incr easing frequency . 
In the case of a point source (i.e. 6 S < Op) (e.g Goodman! 1997 : 
Walkerl ll998l) : 



100 



(2) 



where vo is the transition frequency between weak and strong 
scattering. We fitted this functional form to the three data points 
allowing vo to vary as a free parameter. The resulting fit is 
shown by the dashed line in Fig. [7] The equation describing 
this fit is m = (42.1 ± 3.7) (iip;) 17/12 . 

If the source cannot be regarded as point-like, (i.e. 8 S > Of) 
the modulation index decreas es with increasing source size (e.g 
lGoodman|[l997l: IWalkerlfl998ll : 



m = 100 



(3) 



Since the source size of most extragalactic sources is much 
larger than the Fresnel scale, the latter equation can be used 
to describe the frequency dependence of the modulation in- 
dex of IDV sources. We also fitted this more realistic see- 
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Fig. 7. Frequency dependence of the average modulation in- 
dices. The dashed grey line represents a fit to the data, where 
we assume a source size smaller than or equal to the Fresnel 



scale. Thus m 



-17/12 



The black, solid line shows a fit to 



a more realistic scenario, where we assume that the source is 
extended and that its size changes with frequency following a 
power law dependence. 



nario, where we assumed that 6 S follows a power law depen- 
dence with the frequency. This fit is shown by the solid line in 
Fig. [7] The equation describing this fit is m = (21.5 ± 0.7) • 

/ \-0.77±0.03 ._. 

ItcSHzJ ' Using equation|3] we can calculate the depen- 

/ \-l.l±0.04 

dence of the source size on frequency: Q s ~ I . * m I 
This is in good agreemen t, with the standard jet model (e.g. 
Blandford & Konigl Il979h . Here a constant brightness tem- 
perature is assumed for the self-absorbed jet base: T% ~ 
S /(v 2 f9j) ~ const. Thus, one obtains for a source with flat spec- 
trum a source size, which is inversely proportional to the ob- 
serving frequency. 

Since the modulation index increases with decreasing fre- 
quency, we can give only an upper limit on the transition fre- 
quency between weak and strong scattering of < 2.70 GHz. 
From the annual modulation model we obtained the scintilla- 
tion length scale, so, which is the product of the source size and 
the distance to the scattering screen. With known variability in- 
dex m from the fit and so, we now can calculate the distance of 
the scattering screen and the source size at 4.85 GHz. Since we 
have only an upper limit to the transition frequency, the result- 
ing distance must be regarded as lower a limit, the source size 
as an upper limit. In Table|3J we summarize these values for the 
different annual modulation models. In every case, the source 
sizes are larger than the Fresnel angles, as expected for an ex- 
tended extragalactic source and thus confirming the assumption 
that the scintillation is quenched by a finite source extent. 

The estimates of the source size compare well with the 
core size measured by VLBI (CJF survey, epoch 1994.702, 
0.13 mas). The equipartition source size is about a factor of 
~ 3 larger than the ones obtained here. Doppler factors of 300 
to 600 would be needed to reduced these to the scintillating 



Table 5. Lower limits to the screen distance (Col. 2) and upper 
limits to the source size (Col. 3). In the calculations we used 
the scattering length scale as derived from the annual modula- 
tion model, the frequency dependence of the source size, which 
was derived from the frequency dependence of the modulation 
index, and we assumed 2.70 GHz as the upper limit to the tran- 
sition frequency between weak and strong scattering. In Col. 1, 
the letters refer to the different annual modulation curves, see 
Tableland Fig. [6] In Col. 4, for comparison, we list the Fresnel 
angle at 4.85 GHz at the distance of the screen. Col. 5 gives 
the estimated scattering measure for a transition frequency of 
2.70 GHz. 
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sizes. Another explanation is that the scintillating component 
is only the core, while the spectrum and the turnover frequency 
of J 1 128+5925 are influenced by the non-scintillating compo- 
nent as well. Therefore, the derived equipartition size might not 
represent the scin tillating size. 



According to iTavlor & Cordesl dl993l) . the transition fre- 
quency between weak and strong scattering is vo = 185 - 
SM 6/17 ( T ^) 5/1? , where the Scattering Measure (SM) is de- 
fined by the line of sight integral of the amplitude of the elec- 
tron density fluctuations {C 2 N ) of the turbulent ISM: SM = 

Cffdz, where d is the path length through the scattering 
medium. Thus, knowing the distance, we can also calculate the 
SM. In Table|5]we show the SM derived in this way. 



4.3.1 . Cross-correlation between the frequencies 

In Figure [8] simultaneously measured multi-frequency light- 
curves are displayed. A correlation between the 2.70 GHz and 
4.85 GHz data is obvious. This is consistent with the picture 
that the variations at both frequencies are due to scintillation 
in the same, weak regime. This correlation is quantified by 
the cross-correlation functions shown in Fig. [9] The cross- 
correlation functions reveal a (insignificant) time lag between 
the 2.70 GHz and the 4.85 GHz light-curve of (0.05 + 0.05) day 
in December 2005 and in February 2006. In April 2006, the 
time lag is larger (-0.50 ±0.05) day, showing that the 2.70 GHz 
light-curve is leading. 

The cross-correlation function of the 4.85 GHz and 
10.45 GHz light-curves of May 2005 does not show a signif- 
icant peak (its maximum value is 0.4 ± 0.2). This may be due 
to the sparseness of the data at 10.45 GHz. For the observations 
in December 2005 the correlation is higher, peaking at a time 
lag 0.00 ± 0.05 day with a correlation coefficient of ~ 0.8. 

The differences between the observed time lags may be 
caused by a frequency dependent structure of the scattering 
medium or of the scintillating source. It is possible that the ax- 
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Fig. 8. Light-curves of J 1128+5925 at those dates when measurements were carried out at more than one frequency quasi- 
simultaneously. On the abscissa the observing time in days is displayed, the starting dates are given in the upper right corner of 
each plot. On the ordinate the normalized flux densities are given. Filled grey circles represent the 4.85 GHz measurements, open 
black squares represent the 2.70 GHz measurements, black filled triangles represent 10.45 GHz measurements. 



ial ratio or the position angle (or both) of the anisotropy shows a 
certain frequency dependence and at different times of the year, 
the scintillation patterns are displaced relative to each other. 
According to the anisotropic annual modulation scenario, the 
relative velocity vector points almost at the position angle of 
the major axis of the scintillation pattern in April at 4.85 GHz. 
Therefore the variability was long at this time. In December, 
however, the velocity vector was almost perpendicular to the 
velocity vector in April, pointing closer to the direction of the 
minor axis of the scintillation pattern. The different time lags 
may arise because the ellipticity of the scintillation pattern at 
2.70 GHz is larger than at 4.85 GHz; the major axes at the two 
frequencies differ by more than the minor axes. 



The larger and better defined time lag in April 2006 might 
also reflect frequency dependent size variations in an evolving 
radio source. This would require that the VLBI structure of the 
source changes on a time scale of ~ 0.5 yrs, which is not un- 
reasonable. 



5. Discussion and Summary 

J 1 128+5925 is a new highly variable IDV source, which varies 
on time-scales of 0.2 days to 1.7 days with peak amplitudes of 
up to ~ 30 % and a modulation index in order of 10 %. 

Adopting the usual light-travel time argument and by this a 
source intrinsic interpretations of the variability, the measured 
time-scales lead to brightness temperatures of ~ 10 19 ~ 20 K. 
This requires Doppler factors of at least a couple of hundred 
to reduce the brightness temperature to the inverse Compton or 
to the equipartition limit. 

The variability time-scale of J 1 128+5925 is different at the 
different observing epochs and varies systematically over the 
year. We determine the variability time-scale by three differ- 
ent methods, which ensures the robustness of the effect. We 
explored the possibility if these changes can be described by 
the so-called annual modulation scenario, where variability is 
caused by interstellar scintillation, and where the scintillation 
time-scale is modulated by the orbital motion of the Earth 
around the Sun. We found that such annual modulation fits the 
data, if some anisotropy is included. Similar models were used 
to describe the annual modulation observed in the much more 
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Fig. 9. Cross-correlation functions of the light-curves of 
J 1128+5925 at 2.70 GHz and 4.85 GHz for three different 
epochs. A negative time-lag indicates that the lower frequency 
data are le ading. The cross-correla tion was calculated using the 
method of lEdelson & Krolild dl988l) . 



rapid scintillators J 1819 +3845 dDennett-Thorpe & de Bruynl 
20031) and PKS 1257-326 dBignall et al.ll2006l) . 

Whether the anisotropy originates from the source or from 
the scattering plasma is at present not clear. VLBI observa- 
tion (CJF survey) of J 1128+5925 show a slightly elongated 
core-jet structure, which is oriented at a position angle of 

(91... 100)°. This is in notable agreement with the position 

angle of the major axis of the anisotropy, as derived from the 
fitting of the annual modulation ( 90°). Since the VLBI ob- 
servations were performed more than ten years before our IDV 
measurements, it is unclear if the orientation of the source is 
still the same. Further VLBI monitoring is necessary to search 
for a possible correlation between source structure, its orienta- 
tion and the anisotropy parameters of scintillation models. 

The frequency dependence of the variations (modulation 
index decreasing with increasing frequency) is in agreement 
with the theory of weak ISS. Using the frequency dependence 
of the modulation indices, we derived a source size, which is 
inversely proportional to the frequency. This is in agreement 
with the standard jet model usually applied for a flat-spectrum 
radio sources (AGN). 

Our observations yield an upper limit for the transition fre- 
quency between weak and strong scattering of 2.70 GHz. Using 
the scintillation length scale from the annual modulation mod- 
els and the frequency dependent source size (the size is propor- 
tional to v _1 ), we derived a lower limit to the distance of the 
scattering screen of > 50 pc for an isotropic scintillation model 
and > lOOpc for an anisotropic model. 

VLBI observations of extragalactic radio sources yield 
source sizes which often are sufficiently small for showing in- 
terstellar scintillation. Such sources also should show an annual 
modulation in their variability pattern, if it is caused by scintil- 
lation. However up to now, a seasonal cycle has been unam- 
biguously identified only in two sources, both of which belong 



to the class of extremely fast IDV sources (rapid scin tillators 
Dennett-Thorpe & de Bruynll2003l: iBignall et alll2006h . In the 
case of the slower, classical type II IDV sources, the detection 
of annual modulation is more difficult, since the variability is 
slower and can be less well extracted from time limited data 
trains. Furthermore source-intrinsic structural variations also 
exist in many AGN, on time-scales of months to weeks. This 
can smear out or even d estroy a seasonal cycle. This might 



be the case in 0716+714 dBach et alj|2006l) . Variations of the 



intrinsic source structure or of the scattering medium could 
also explain temporal variations in the IDV pattern, and per- 
haps even a sudden stop of such variability. Such an effect may 
be seen in B 0917+6 24, for which an annual modulation was 
previ ously proposed, dRickett et al hooit iJauncev & Macquarn 
20011). but later could not be confirmed due t o the ceasing of 
IDV dKraus et al.ll 19991: iFuhrmann et"ai]|2002h . In fact a VLBI 
monitoring performed in parallel, shows structural variations in 
the core region of B 0917+624, which may be related to the ob- 
served decrease of the IDV amplitude and prolongation of the 
variab ility time scale dKrichbaum et al. 2002; Bernh art et al 
2006b . We n ote that episodic IDV wa s also observed in PKS 



B 0405-385 dKedziora-Chudczerll2006T) and similarly was in- 
terpreted as being either due to variations in the physical prop- 
erties of the scattering screen or a morphological evolution of 
the source. 

To investigate whether source intrinsic variations play a 
role in the variability of J 1128+5925 further VLBI observa- 
tions have to be carried out. During the 1 .5 years of observa- 
tions, the source experienced a 20% increase followed by a 
similar decrease in its flux density at 4.85 GHz. This could in- 
dicate variations of the intrinsic source structure. Future kine- 
matic VLBI studies of J 1128+5925 can also help to derive a 
more accurate estimate of the Doppler factor, and by this of the 
source size and the resulting "quenched" scintillation. Together 
with the scintillation length-scale derived from the annual mod- 
ulation model, this would give more stringent limits to the dis- 
tance of to the scattering screen. 

In principle one needs to observe the source regularly for a 
longer period of at least two or three years to unambiguously 
identify an annual cycle in the scintillation. The observations 
of J 1128+5925 presented here cover just 1.5 years. Only two 
observations were carried out at the same time of year, with 
one year separation (2004 December and 2005 December). 
The variability time-scales extracted from these observations 
are in very good agreement with the predictions of an annual 
modulation model. A prolongation of the variability time-scale 
can also be caused by an increase of the size of the scintil- 
lating component in the source, for example because of the 
ejection of a new jet component. Therefore, observation of the 
source during the time when we expect slower variability (here: 
in April-May and in August-September) is highly desirable. 
Over-subscription and logistical constraints at large radio tele- 
scopes make such observations difficult. It is therefore useful 
(and necessary) to share the observational load between differ- 
ent telescopes. The combination of multi-frequency observa- 
tions at two complementary telescopes (Effelsberg & Urumqi) 
provides an efficient way to observe and detect systematic vari- 



K. E. Gabanyi et al.: J 1 128+5925: Annual Modulation in a new IDV source 



15 



ations in the intra-day variability of radio sources, which vary 
slower than the so called rapid scintillators. 
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